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Abstract 

In in this paper we show how using D.A. it is found a simple change of variables (c.v.) that brings 
us to obtain differential equations simpler than the original one. In a pedagogical way (at least we 
try to do that) and in order to make see that each c.v. corresponds to an invariant solution (induced 
by a symmetry) or a particular solution, we compare (with all the tedious details, i.e. calculations) 
the proposed method with the Lie method. The method is checked even in odes that do not admit 
symmetries. 
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1 Introduction 

D.A. has usually been employed in different areas (fields) such as engineering problems, fluid mechanics 
etc... and these problems are always described by partial differential equations (pde) (see [1]-[10]). This 
method ("tactic") helps us to reduce the number of quantities that appear into an equation and to 
obtain ordinary differential equations (ode). We would like to point out (emphasize) that this tool is 
more effective if one practices the spatial discrimination, such tactic allows us to obtain better results 
than with the standard application of D.A. (see [11]). Knowing that D.A. works well in pde we would 
like to extend this method to the study of ode (the first order in this case) in a systematic way. There 
are in the literature previous work in this direction for ode of first order (see [12]-[13] . 

The method of the Lie groups it has showed as a very useful tool in order to solve nonlinear equations 
(nl-ode) as well as pde (see [14]-[27] ). Nevertheless, when one is studying ode of first order its application 
results very complicated (very tedious) if one has not a computer algebra package since if one decides 
to look for the possible symmetries of an ode with pencil and paper this task may be turned very 
exhausting. However, if we know that an ode admits a concrete symmetry, then it is a trivial issue to 
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find new variables which allows us to rewrite the ode in quadratures or as we will see in this paper to 
obtain an ode with separating variables. 

Our purpose in this work is to explain, through examples, how D.A. works in order to find these 
changes of variables (c.v.) in a trivial way, i.e. without the knowledge of the symmetries of the ode 
under study. The idea is as follows. When we are studying an ode form the dimensional point of view, 
we must require that such ode verifies the principle of dimensional homogeneity (pdh) i.e. that each 
term within the equation have the same dimensions, for example, speed, or energy density. To clarify 

y 

this concept we consider the following ode, y' = — \-x, where each term must have dimensions of y' i.e. 

x 

W] = V X ~ X i where [•] stands for dimensional equation of the quantity •. As we can see, this ode does 
not verify the pdh, since the term x has dimensions of x, i.e. [x] = x. In order to do that this equation 
verifies the pdh we need to introduce dimensional constants, in such a way that after rewrite the ode 
with these constants the ode now verifies such principle i.e. in this example and as we can see easily 

if we consider only one constant a, such that [a] = yx~ 2 , we make that the ode y' = — h ax, verifies 

x 

the pdh. We would like to emphasize that this situation does not appear (arise) when one is studying 
physical or engineering problems since (as it is supposed) that such problems (equations) verify the 
pdh and we do not need to introduce new dimensional constants, that must have physical meaning, for 
example, the viscosity coefficient, etc.... 

Precisely this dimensional constant suggests us the c.v. (x, y{x)) i — ► (t, u(t)) where (t = x,u = yx~ 2 ) 
in such a way that rewriting the original ode in these new variables it is obtained a new ode with sepa- 
rating variables: u + u't = 1. The reason is the following. We know from the Lie group theory that if 
it is known a symmetry of an ode then this symmetry bring us through a c.v. to obtain a simpler ode 
(a quadrature or a ode with separating variables) and therefore the solution is found in a closed form. 
To find this c.v. it is used the invariants that generate each symmetry, i.e. the first principle is that 
it is useful to pass to new coordinates such that one of the coordinate functions is an invariant of the 
group. After such transformation it often (but not always) happens that the variables separate and the 
equation can be solved in closed form. We must stress that taking the new dependent variables to be 
an invariant of the group does not guarantee the separation of variables. The choice of the independent 
variable is also very important. 

One must note that D.A. gives us this invariant (or at least a particular solution). It is observed in 
our example that the solution y = x 2 , (suggested by D.A.) is a particular solution for this ode, but if we 
study this ode form the Lie theory point of view it is obtained that such ode has the scaling symmetry 
X = xd x + 2yd y and therefore the solution y = x 2 , is furthermore an invariant solution (generated by 
X), for this reason such c.v. works well. As we will show in the paper, the c.v. that D.A. induces is 
the same than the generated one by the symmetries of the ode. In order to make see this fact, we will 
solve each example by D.A. as well as by the Lie group technique (abusing of the trivial calculations), 
calculating the symmetries and their corresponding invariants and c.v.. 

However D.A. has a little (or rather big) drawback, it only gives us relationship of type power i.e. 
y = x n , with n £ R, and we cannot obtain relationship as y = (x + l) n . Nevertheless, and as we will 
see in the examples, although D.A. does not give us an invariant solution it will be sufficient that it 
provides us a particular solution in order to obtain a c.v. which brings us to a simpler ode than the 
original one. This task will be very useful in the case of studying Riccati and Abel odes. Furthermore, 
this "tactic" is valid even in the case in which the ode has no symmetries but unfortunately one always 
founds examples in which any tactic does not work. 

The paper is organized as follows: In section 2 we describe (with all the tedious and superfluous 
calculations), through two examples, how the proposed method works by comparing it with the Lie 
method in order to show which are the symmetries and their corresponding invariants. In section 3 we 
will show several examples beginning with a Bernoulli ode and following with two Riccati ode and three 
Abel ode (see [28]-[32]). We think that these kind of ode are the most difficult to solve and therefore 
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any simplification is welcome. In section 4 we show two Abel odes that do not admit symmetries (see 
[33]-[37]). In the first of them we use the D.A. to obtain a particular solution that allows us to obtain a 
Bernoulli ode but in the second one we are not able to find any particular solution and therefore a c.v. 
which brings us to obtain a simpler ode. We end with some conclusion as well as pointing out some of 
the limitations of the proposed tactic. 

We would like to emphasize the pedagogical character of this paper (at least we try) for this reason we 
have abused of superfluous calculations and we have omitted some technical details since the supposed 
audience are mainly engineers and/or physicists but not for mathematicians. 

2 The method 

In this section we will explain how the D.A. works in order to solve odes through two simple examples. 
The main idea is to introduce dimensional constants that make the ode under study verify the principle 
of dimensional homogeneity. These dimensional constants help us to find c.v. which bring us to obtain 
ode with separating variables. These c.v. obtained through D.A. correspond to invariant solution (or at 
least to a particular solution) and therefore they are induced by one symmetry. We compare our tactic 
with the Lie one. 



Example 1 Solve the homogeneous equation 

(x 2 + y 2 ) dx = 2xydy. 



(1) 



Solution. In order to solve it we will use three different ways, the traditional, the dimensional 
and the Lie method. 

Traditional method. Making the c.v. u = y/x we have: 



1 - u 1 



u 



2udu dx 

x 



ln(l — u 2 ) = lnx ==>• y 2 = x 2 + x, 



2ux 1 — u 2 

Dimensional Analysis. We go next to consider eq. (1), written as follows 

/ _ (a 2 x 2 + y 2 ) 



V 



2xy 



(2) 



(3) 



where the dimensional constant a, makes the ode verify the dimensional principle of homogeneity (d.p.h.) 
if 



x y, 



(4) 



where [•] stands for the dimensional equation of the quantity •. 

Applying the Pi theorem we obtain the dimensionless variables that help us to simplify the original 
ode. Therefore taking into account the following dimensional matrix we take 





y 


X 


a 


X 





1 


-1 


y 


i 





1 



ax 

y 



y = ax, 



(5) 



as we will see later this solution is a particular solution of this ode. We would like to emphasize that as 
we have only needed one constant then the equation is scale-invariant in such a way that the generator 
of this group is X = xd x + yd y as it is observed from eq. (5). This fact will be probed studying this 
equation under the Lie group tactic, see below. 
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In this way the new variables are (t,u(t)) : 



t = x, u(t) = a^-j , =>- (^x = t, y= a-^y ) 



this change of variables brings us to rewrite eq. (3) as follows: 

u' 1 u 1 



u(l-u 2 ) 2t vT^ 5 

and hence 



Int, (7) 



y 2 = a 2 x 2 + dec. (8) 

As we can see in this trivial example, the D.A. induces a c.v. which helps us to obtain an ode simpler 
than the original one. 

We can also think in the following way 

ay ' = b Yy +a b (9) 



where 



and hence 



[a]=x, [b]=y z x- 1 , (10) 
t = — , u(t) = j— ) (^x = at, y = \J abtu{t)^ , (11) 



therefore eq. (9) yields 

in this way we obtain the solution 



v! = i=> u = t + C 1 , (12) 



y 2 = -x 2 + ClX, (13) 

a 



once we have obtained the solution, the constants a, b are makig equal to 1, i.e. a = b = 1. 
In the same way we may consider the following change of variables 



u(t) = -j= j =>• (x = at, y = u(t)Vabt) , (14) 



x , y 

u{t) = 
a y/bx 

which brings us to the following ode 

2u'u = 1 u 2 - t + d = 0, (15) 

and therefore we obtain again the solution (13). As we will see below all these c.v. are generated by 
their corresponding (respective) symmetry. 

As we can see, this last change of variables is better than the first and the tactic is the same: to 
introduce dimensional constants that make the equation dimensional homogeneous. 

Lie Method. In order to find the symmetry generator of a ode 

y' = f(x,y), (16) 

we need to solve the following pde 

Vx + (Vy ~ £x)f ~ iyf ~ y)f x - rj(x, y)f y = 0, (17) 
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where 

/. = f, /. = f (IB) 
ax ay 

In this case we have to solve 



2 xy 4 x z y z y 2x z y x 2xy z 



we find that [X = £d x + rjd y ): 



x I x^ 1 y 

Xi = — d y , X 2 = ( -— + g ) X 3 = xd x + ydy, 

^4 = d x + ^-d y , X 5 = (x 2 + y 2 )d x + 2xyd y , (20) 

observing that X§ is a trivial symmetry. 

Each symmetry induces a change of variable (canonical variables) which are obtained through the 
following formula 

Xt = 0, Xu = 1. (21) 
In this way the change of variables that induces the field X^ = d x + ^d y is: 

t = ^,«(t) = x) => (x = u(t),y = tyftfify . (22) 
in such a way that eq. (1) yields: 

it' = 2t =► it(t) = i 2 + Ci , (23) 

and therefore 

y 2 

x = — + Ci, (24) 

x 

which is the same solution. The c.v. induced by this symmetry is similar to the obtained one in (14). 
The c.v. that induces the field X\ = ^jd y is the following one 



t = x, u(t) = \ => (x = t, y = y/u(t)tj , 



(25) 



therefore eq.(l) yields 

u' = b 2 , (26) 

this c.v. is the same than the obtained one in (11). 

For example the c.v. that induces the field X3 = xd x + yd y is 

t = | , u(t ) = In x) => (x = e u W , y = te u W ) , (27) 

which brings us to the following ode : 

2t 

u = =>u(t) = -ln(-l + t) -ln(l + t) + Ci, (28) 

r z — 1 



now writing the solution in the original variables we take 



= Cl + 1 °l, fa-x)fa + x) ,)- (29) 
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Now if we consider the invariants that induce each symmetry 



dx dy , dy 77 . . 

-r = — , 1 — ► y ■= -T- = 7, (30) 



we take that: 



X 1 ^I 1 = x, X 2 ^I 2 = x, X 3 ^/ 3 = -, X 4 ^/ 4 = ^=. (31) 



X x X 



For example the symmetry X3 = xd x + yd y generates the following invariant: 
dx dy dx dy y 

—— = — ==>• — = — ==>• In x = my ==>• 13 = — 1 — ► y = ax, (62) 
5 r\ x y x 

this would be the solution that suggests us precisely the direct use of the Pi theorem (if the ode is scale 
invariant, as in this case, the solution obtained applying by the Pi theorem coincides with the invariant 
solution that generates the scale symmetry, in this case X3). We see that we only obtain a particular 
solution, but that this is invariant, in fact if we think about the ode as a dynamical system we see that 
the fixed point of such equation would be precisely the solution y = ±x. ■ 



Example 2 Solve the linear ode 

(l-x 2 )y' + xy = 1. (33) 

Solution. Solution through D.A. Our first step will be to introduce dimensional constants in 
such a way that eq. (33) verifies the principle of dimensional homogeneity. In this way we rewrite the 
equation as follows 

{A - x 2 ) y' + xy = B, (34) 

where 

[A] = X 2 , [B] = XY, (35) 

i.e. 

x 2 xy 

711 = ~A ^ = ~B' ^ ' 

where 

* = iMt)' (37) 

being (p a unknown function. It is possible to make the following assumption 

V = Wx {^) , n€R. (38) 

Since y = x _1 is not a particular solution of (33) we need to combine the 7r — monomias in order to 
find a particular solution, the simplest way is as follows: 

t = ^, u (*) = ^|) ( x = "/Ali y = ju(t)^A^ (39) 
observing that 

U (t) =TT 2 -^\ (40) 

where y = cx, derived from tt 2 ■ vr^ -1 is a particular solution of (33). 
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Remark 2.1 (Recipe) When we have two ir—monomia, we have to check if they induce any particular 
(or invariant) solution. If they are not particular solutions then we must combine them in order to obtain 
such solution, this solution usually is obtained by combining them in a very simple way. 

In this way our ode is written now as follows: 

' u = 1 + Ci — - — , (41) 



u-1 t(t-l) t 
and in the original variables the solution of eq. (33) yields: 

y = ^(x±CWA-x^, (42) 

where C\ is an integration constant. Now making A = B = 1, we have the ordinary solution to eq. 
(33). _ 

Lie method. Following the standard procedure we have to solve the pde: 

(t)v — £x) (1 — xy) £ v (l — xy) 2 .. y 2(l — xy)x. ri(x,y)x . , , 

rix + K ' V p V 2 --^ 2§- V) -i-^T + -)i ^~ + 2 = °» 43 

1 — (1 — ar J z 1 — x z (1 — x z ) z 1 — x z 



which solutions are: 



X 1 = V^Tx~ 2 d y , X 2 = (-y + x)d y , X3=( i2X \j! 1) ^ i + y i) +X) )dy> (44) 



and their corresponding invariants are: 

X 1 ^I 1 =x, X 2 ^I 2 = x, X 3 ^I 3 = x. (45) 
For example the symmetry X\ generates the following c.v.: 



V-l + x< 



t = x, u(t) = )=>(x = t, y = u(t)V-l + t 2 ) , (46) 



therefore eq. (33) is written as: 



u> = - (47) 



which solution is: 



{ -l + t 2 f 2 (-l+t 2 ) 3/2 
4(t + l) 2 4 (* - 1) 



u(t) = - K — / +^— + (48) 



and hence 

y C x (x 2 - 1) +xVx 2 - 1 



y = x ± cVx - lVx + 1. (49) 



Vx^T x 2 -l 
Now, if for example we consider the symmetry X 2 then we have: 

/ te u ® — l\ 
(t = x, u(t) = - ln(-y + x)) I x = t, y = ^ — I, (50) 



therefore: 
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finding in this way that its solution is: 

u(i) = -iln(i-l)-^ln(t + l) + C 1 , (52) 
hence in the original variables the solution yields: 

- ln(-y + x) = -- ln(x - 1) - - ln(x + 1) + C u (53) 
which after a simple simplification it yields: 

y = x ± cVaJ — + 1, (54) 

as we already know. ■ 



With these two simple examples we have tried to show how the D.A. works in order to introduce 
c.v. which brings us to obtain simpler odes than the original one. As we have emphasized in the recipe, 
the trick is to look for a particular solution. Sometimes this particular solution will be furthermore 
invariant solution (induced by a concrete symmetry). If this is the case, then our ode will be reduced 
to an ode with separating variables but if this solution is only a particular then, as we will see below, 
we have not any guarantee of reducing our ode to an ode with separate variables, nevertheless we will 
obtain a simpler ode than the original one. In the next section we will study some examples. 

3 Examples 

In this section we will apply our pedestrian method to different odes, beginning with a Bernoulli ode. 
We go next to employ the method to solve Abel odes as well as Riccati odes since as anyone knows 
these kind of odes are truly very difficult. 



Example 3 Solve the following Bernoulli ode. 

, y (x + 1) y 2 
V= -—l — ' (55) 

Remark 3.1 Historically it was Bernoulli the first person who introduced c.v. in order to solve odes 
(now bearing his name). He managed to reduce this equation to a simpler linear equation. 

Solution. It is observed that with 

[a] = X, [b]=X- 2 Y-\ (56) 

eq. (55) yields d.h. i.e. 

x + a 2 
in this way it is obtained the following variables 

t = ~a> u{t) = ^Ll)^{ X = at > y = JbTu)' (58) 
with this change of variables eq. (55) yields 

2u + 2t{l + t)u = {t + lf , (59) 
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which is linear and its solution is: 



and therefore in the original variables the solution yields 

V= (bx + C^ix + a)' (61) 

in this case D.A. does not bring us to obtain a good change of variables but helps us to obtain a simpler 
ode than the original one. We can try now with the following c.v. 

t = x, u(t) = ±- J =► (* = t, y=^), (62) 

which brings us to obtain a linear ode 

t(t + 1) (2u't + 4u) = 2t 2 u + t 2 + 2t + 1, (63) 

but with this tactic we are not able to obtain a simpler ode. 

As we can see D.A. induces us to the c.v. but y = — or y = — , are not particular solutions of 

(55). Nevertheless solutions as y = or y = ^ are particular and invariant solutions for this 

x + i {x + \y 

reason these c.v. induce us to obtain an ode in separate variables, but unfortunately D.A. is unable to 
construct such c.v.. For this reason we can only reduced our ode to a linear ode as in the theoretical 
case i.e. employing the theoretical method purposed by Bernoulli. 

If we study this ode under the Lie group method, it is observed that eq. (55) admits the following 
symmetries obtained from: 



y {x + l)y 2 \ t ( y (x + l)y 2 \ 2 



Vx + (rj y - -—— - ty 



x + l 2 J x + l 2 



y i 2 \ / i 

( X ^f-2 y ) - 7] {- y - Xy -- + l 



< ( ; ^ 2 " *y 2 ) -v\-y-xy- z~r ) = ( 64 ) 



y 2 (x + l) y (l + 2x)y 
Xi = d y , X 2 = [x + 1) d x - 2yd y , X 3 = d x - + ^ d y , X A = xd x u~^f\ d y 

(65) 

and which respective invariants are: 

h = x, I 2 =y(x + l) 2 , h = y(x + l), h = yx(x + l), 
that induces the following canonical variables. For example from X\ : 

t = x , u(t) = j— — -)=^(x = t, y= — — -V (66) 

y(x + a)J \ u(t + a) ) 

in such a way that eq. (55) yields 

u' = -b=> y = bt + Ci, (67) 
and in the original variables i.e. (x, y) this solution yields 



V (bx + Ci)(x + a)' (68) 



as we already know. 
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Example 4 Solve the following Riccati ode 

i 2 2y 1 . , 

y =xy 7. (69) 

Solution. It is observed that introducing the following dimensional constants eq. (69) verifies the 
principle of dimensional homogeneity (p.d.h.), 

, xy 2 2y a ... 

y' = — - -3, 70 

where [a] = X 2 Y. As we can see if we need only one constant then our equation is scale invariant and 
therefore this constant indicates us that the generator of this symmetry is X = —2xd x + yd y as we will 
see below. Therefore we have the following variables 

t = x , u(t) = —2 J => (x = t, y = J , (71) 
yx z J \ tu z / 

with this new variables eq. (69) yields 

u' 1 



In t + arctanh(u) + d = 0, (72) 



n 2 - 1 i 

and in the original variables this solutions is written as: 

lnx + arctanh(-^r) + C\ = 0. (73) 
yx z 

We would like to emphasize that in this case, as we have only one dimensional constant, we can find 
the particular solution y = ax~ 2 applying the Pi theorem, and that such solution (invariant particular 
solution) is induced by the scale symmetry. 

Eq. (69) equation has the following symmetries: 



and hence 



12 1 

Xx = ~d x - -^dy, X 2 = x 4 (y + ^fdy, , ^3 = xd x - 2y d y 

X A = (-^ + y 2 x 2 ^jd y , Xs = x 3 d x - (2 + 4yx 2 )d y . 



(75) 



and which corresponding invariants are: 



h = yX x2 1 1 I 2 = h = x, I 3 = yx 2 , I 5 = yx 2 (x 2 + l) . (76) 



The method of canonical variables brings us to obtain the following ode, for example, for the sym- 
metry X5, we have 

t = x 2 {yx 2 + l), U (t) = -J_), (77) 

2 x z J 

ut = ^ =>u(t) = ~\ + Ci, (78) 
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and hence 

2x 2 x 2 (p 2 + l)^ 11 ^ 

as we already know. 

The transformation induced by X3 (scaling symmetry) is the following one: 

(t = yx 2 , u(t)=hx(x))=>(x = e u ®, y = 7^2) > ( 80 ) 

ut = -= => ti(£) = — arctanh(i) + C\, (81) 

— 1 

therefore 

ln(x) = — arctanh(y x 2 ) + C\, (82) 
which looks a little different than the other solution. ■ 



Example 5 Solve the following Riccati ode: 

y = H -3- (83) 

Solution. The ode verifies the p.d.h. if we introduce the following dimensional constant: 

/ V + a y 2 

y = + &-,, (84) 

x x A 

where [a] = Y, [b] = X 2 Y~ 1 . In this way it is obtained the c.v. 

t = , u(t) = — ^ ==$■ (x = V abtu 2 , y = au) , (85) 
and therefore eq. (83) yields 

u 

u + = 

(ut + t + 1 



u ut , . 

-,+t= 7-1 r~TT~iT ' (86) 



which is a Bernoulli ode and its solution is: 

1 



y/i (^arctan y/i + C x J =0, (87) 
hence in the original variables the solution is 



a 

y 



0. 



It is observed that we have the particular solution y = Cx, obtained from the relationship 
Since this c.v. is not very good for our purposes we may now proceed as follows: 

x 

where 



U2V_ 
,-2 



axy r = ay + 1 + b 
a]=y-\ [b}=y~ 1 x (90) 
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in such a way that these constants induce the following change of variables: 

a x \ ( b t 



which brings us to obtain 



- + arctan(u) + d = (92) 



u 2 + 1 t 2 t 
and hence 

— - + arctan(^) + C\ = (93) 
ax by 

which is a better approximation. In the same way we can obtain a quadrature simply considering t = — 

ax 

and u{t) = j^-. in this way it is obtained u' = u 2 + 1. 
The Lie method bring us to solve the following pde 

Vx + (Vy ~ £x)( ~ + -3) " tyi— + -3) " £( -2- -—4-)- V(~ + -3 ) = 0, (94) 

and which solutions are: 



sin(| ) 2 (ysin(^) - co s(|) a;) 2 
rr(-l+cos(i) 



^3 = I ' , „ dy, X 4 = x 2 d x + xyd y (95) 



and which invariants are: 



h = h=h = x, h = -. (96) 



The symmetry Xi = ( x + — J 9 y generates the following c.v.: 



(t = x, u{t) = arctan y~Jj =^ (y = tan(u(t))i, x = t) , (97) 

u t = p^u{t) = -- + C 1 , (98) 

and hence the solution is 

arctan (^) = h C h (99) 



ay x 
as we already know. 

We would like to emphasize that this ode is not scale invariant and nevertheless we have been able 
to obtain a c.v. that reduces the ode to a quadrature. ■ 



Example 6 Solve the following Abel ode 



1-v 2 

y' = w — + 1 (100) 

xy 
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Solution. Introducing the following dimensional constants, we make that eq. (100) verifies the 
p.d.h. 

(xy) y' = -y 2 + bxy + a, (101) 
where [b] = [a] = y 2 . Therefore the c.v. that suggests the D.A. is the following one: 

t = u (t) = y-\^( x = -fL, y = ^) , (102) 

a a I V bwua I 



hence eq. (100) yields 

u (2 + t-u)t = 2u(l + t-u), (103) 



and its solution is: 



U= -2t + 21n(l + 0-Ci ' (104) 



in the original variables it yields: 

b 2 x 2 bxy bxy 



— £+ln(l + — + (105) 
2a a a 

Once again we emphasize that the particular solution, y = Cx" 1 , (in this case invariant solution , 
see below) has been obtained from the relationship 

Alternatively we can try the following c.v. Writing the original Abel ode, eq.(100) in the following 
form: 

y 2 

axyy' = h xy + b, (106) 

a 



where [a] = XY , [b] = XY, therefore we find the next c.v. 



t = ^, u(t) = ( x = v 7 ^, y = -2=) 

b abj V Vabu' 



that brings us to rewrite eq. (106) as follows: 



2t 

u=- =^-u = 2t-21n(l + *) + Ci, (108) 

1 ~\~ L 



and hence in the original variables (x,y) we get: 



4 = 2^-211,(1 + ^) + ^, (109) 
ab b b 

as we already know. We can check that the following c.v. also works well 

x 2 , . xy\ ( i—— bu \ , 

t = — , ut =JL)=>[x = Vabt, y = -= , 110 
ab b ) \ y/abtj 

since in these variables eq. (106) is written as 

2u'u = l + u^t-2u + 2\n{l + u) +Ci = 0, (111) 

and therefore 

^!_2^ + 21n(l + ^) + C 1 = 0, (112) 
ab b b 

obtaining the solution. 
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Applying the Lie method we need to solve the following pde: 

% + + !)-^v( hl ) + 2 + ^- + — =0, (113 

xy x z y x xy z 

which solution is: 

X x = ~d x + ^dy => /i = xy. (114) 
This symmetry induces the following c.v.: 

x 2 \ 

u{t) = - — , t = xy\, (115) 

and therefore in this variables eq. (106) is written as follows: 

u' = — =>u(t) = -t + a1n(a + t) + Ci, (116) 

a + i 

hence in the original variables we get: 

x 2 

— — = — xy + a ln(a + xy) + C\, (117) 

This is another example of an ode that it is not scale invariant, and nevertheless, we have been able 
to reduce it to a quadrature. ■ 



Example 7 Solve the following Abel ode. 

x(y + A)y' = y 2 + 2y + 2x. (118) 

Solution. Following our pedestrian method, we begin by introducing dimensional constants and 
rewriting the odes as follows 

x (y + 4a) y' = y 2 + 2ay + 2bx, (119) 

where [a] = y, and [b] = x _1 y 2 , therefore D.A. suggest us the following c.v. (note that y = y/x, is a 
particular solution of (118)) 

t = t U ^ = i)^{ x = ^ y = ua )> (120) 

in such a way that eq. (118) yields 

t (t + 4) v! = t(u + 2) + 2u, (121) 
which is a linear ode and its solution is: 



u = t + Ciy/t (t + 4), (122) 

hence 

l = f + cj£(£ + 4), (123) 
a bx y bx \bx J 

and simplifying, it is obtained the solution: 

bx 



— = y + C lV / (y 2 + ibx). (124) 
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2 1 

Another particular solution could be found from the following relationship (f) i.e. y = Cx. 
But this particular solution bring us to the following c.v. 

(«=s? »=«-)• < 125) 

which transforms eq. (118) into a Bernoulli ode, 

t 2 + 4) v! = (ut + 2t + 2) (u'i - -u) , (126) 

but this situation is not desirable since it is always more difficult to solve a Bernoulli ode than a linear 
ode. 

Applying the Lie method, following the standard procedure, we need to solve the following pde 

, (Vy-Z x )(y 2 + 2y + 2x) j y (y 2 + 2y + 2x) 2 



x (y + 4) x 2 (y + 4) 2 

2 y 2 + 2y + 2x ( 2y + 2 y 2 + 2y + 2x 

~ a x(y + A)~ x 2 (y + 4) )_7?( x(y + 4)~ x (y + 4) 2 } ' ( } 

and which solutions are: 

f (2y + 4-x)(y-x)\ ( {y 2 + Ax) (y - x) \ , ~ 

Xl = ^ ^4 J a f ^ 2 = ^ x(y + 4) J ^3 = (4x + x ) 8 x +x (y + 4) d y 

(128) 



where their correspond invariants are: 



v + 4 

/i = J 2 = x, J 3 = y —— (129) 
x + 4 



For example the c.v. that induces X% is: 



hence: 



finding that: 



which solution is: 



4(t_l + e (4«(t))) 4e (4u(t)) \ 

y ~ _l + e (4u(t)) ' X ~ _l + e (4u(t)) l> ( 131 ^ 



U ' = 2(2t 2 -3t + iy < 132 > 
u(t) = - ln(i - 1) - - ln(2 t-l) + C u (133) 



therefore, in the original variables it yields: 



\ ln(x) - t ln(4 + x) = - ln(^±i - 1) - \ i n ( 2 ^ + 4 ) _ i) + C ( 134 ) 
4 v ; 4 v ' 2 v 4 + x ; 4 v 4 + x ; v ' 



as we already know. 



Example 8 Solve the Abel ode. 

y > = Cx 3 y 3 + Bxy 2 -A-, (135) 
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Solution. If we rewrite eq. (135) introducing the following dimensional constants, 

y ' = a 2 Cx 3 y 3 + aBxy 2 - A-, (136) 



x 



where A,B,C G K, and [a] = X^Y' 1 . As we can see this ode is scale invariant since we have needed 
introduce only one constant. D.A. suggests us the following c.v. 

(t = x, u(t) = ax 2 y) (x = t, y = -^j , (137) 

in such a way that eq. (135) is written in the following form: 

tu' = u(u 2 + u+l) , (138) 

and its solution is: 

In t + i In (u 2 + u + 1) + ^ arctan ( ( + ^\ V^j - kiu + C\ = 0, 



(139) 



therefore in the original variables it yields: 

lnx + ^ln ([ax 2 yf + ax 2 y + l) + arctan (ax 2 y) + ^ - In (ax 2 i/) + Ci = 0. (140) 

In second place, we study eq. (135) 

y ' = Cx 3 y 3 + Bxy 2 - A-, (141) 

x 

with respect to the dimensional base B = {T} . This ode verifies the principle of dimensional homo- 
geneity with respect to this dimensional base. Note that [y] = [jp] = T 2 , and [x] = [H] = T _1 hence 
[y'\ = T 3 . Therefore rewriting the equation in a dimensionless way we find that y oc x~ 2 

But if we study this equation with respect to the dimensional base B = {X, Y} , we need to introduce 
new dimensional constants that make the equation verify the principle of dimensional homogeneity 

y' = aCx 3 y 3 + (3Bxy 2 - A- (142) 

x 

where [a 1 / 2 ] = [(3] = X^Y' 1 , hence 

>y a ^, (143) 

As we can see we have obtained the same solution than in the case of the invariant solution. This is 
because the invariant solution that induces a scaling symmetry is the same as the obtained one through 
the Pi theorem. 

This ode admits the following symmetry (scale-invariant) 

X = xd x -2yd y , => I = x 2 y (144) 
which is a scaling symmetry and it induces the following change of variables, 

r = x 2 y, s(r) =ki(x), x = e s(r) , y = -^pj, (145) 





y 


P 


X 


X 





-2 


1 


Y 


l 


-1 
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which brings us to obtain the next ode in quadratures 



S ' r (Cr 2 + Br + 2- A)' (146) 

and which solution is: 

Inr l ln{Cr 2 + Br + 2-A) garctaph 
S(r) ""^ + 2 (A-2)VB 2 + 4C(^-2) ' 

and hence in the original variables (x,y): 

\n(x 2 y) lln(Cx 4 y 2 + Bx 2 y + 2-A) Baxctimh \y/B*+4C(A-2) y 

Inx = r J l +Ci, (148) 

A- 2 2 A- 2 [A - 2) ^B 2 + 4C{A - 2) 

which is the most general solution for this ode. ■ 



4 Pathological cases. 

In this section we will present two examples of odes that do not admit symmetries (Lie point symmetries) . 
Nevertheless in the first of them D.A. helps us to obtain a simple c.v. that will bring us to obtain 
a simpler ode (through a particular solution). In the second case we will show that unfortunately 
sometimes one finds odes that at this time have no solution, or at least we do not know how to solve 
them. 



Example 9 Solve the Abel ode 

(x 2 y + x 5 -x)y' = xy 2 -(x 4 + l)y. (149) 
Solution. Following our pedestrian method we beginn by introducing dimensional constants 

ax 2 yy' + bx 5 y' — cxy' = dxy 2 — ex 4 y + fy (150) 

where 



[a] = [d] = x~ 2 y-\ [b] = [e] = x~ 5 , [c] = [/] = x~\ (151) 

therefore [c 5 ] = [/ 5 ] = [b] = [e] , having only two dimensional constants. Since y = x~ 2 is not a 
particular solution of eq.(149) then we look for a combination between the monomias finding in this 
way that 

1 ax y a ko\ 

7Tl = — , vr 2 = (152) 

cx c 

where y = x" 1 is a particular solution of (149). The c.v. that induces the D.A. is the following one: 

(t = xy, u(t) = x) =^ (x = u, y = lj , (153) 

hence using these new variables eq. (149) is written now as: 

2u't(t + l) + u(l-t-u 4 ) =0, (154) 
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which is a Bernoulli ode and its solution is 



>tJ-2t -21n(t- 1) + d 
m= v , - (155) 

• v /-2t-21n(t- 1) + Ci 



undoing the c.v. we find that the solution to eq. (149) is: 



' xy\l 1 —2xy — 2 In (xy — 1) + C\ 
x = v (156) 
^-2xy-21n (xy - 1) + C[ 



Applying the Lie method we see that the pde to solve is: 

.2 



y (xy — x 4 — l) y 2 (xy — x 4 — lV 



Vx + (% Sx) / , A i\ S2/ 

x (xy + x 4 - 1) j; ( xy + x 4 - 1) 

y (y — 4x 3 ) y (xy — x 4 — l) y (xy — x 4 — l) (y + 4x 3 ) 



2 



yx(xy + x 4 -l) x 2 (xy + x 4 -l) x ( xy + X A - 1) 

^x(xy + x 4 -l) (xy + x 4 -l) (xy + x*-if) 

but in this case we have not found any solution, i.e. eq. (149) does not admit symmetries. Nevertheless, 
one always may try to find, as if by magic, any c.v. that brings to find a simpler ode. ■ 



In this example we would like to show that there are some odes which are very intractable. Our 
pedestrian method does not work in this case, we have not been able to find any particular solution. 
The Lie method does not work, i.e. this ode does not admit any symmetry and the theoretical methods 
do not work either. At this time we do not know how to solve it. 



Example 10 Try to solve the Abel ode 



y' = -Ax 2 y 3 - By 2 - -, (158) 

x 



where A and B are dimensional constant. 



Solution. In this occasion we already have the dimensional constants 

[A] = x~ 3 y~ 2 , [B] = IT 1 *" 1 , (159) 

in such a way that all the terms of the equation has dimensions of | . We check if one of the dimensional 
relationship induces any particular solution finding that this is not the case. Therefore we go next 
to look for any trivial combination between them but we are not able to find any particular solution. 
Nevertheless we try to obtain any result with the following c.v.: 

/ B 2 1 \ / B 2 At \ , . 

\ t= A? < t) = W X )^\ X= M' y =BW)) (160) 

which bring us to obtain the next ode: 

t 2 u'u + tu + 1 = 0, (161) 
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but unfortunately we have not advanced. 
Another try is the following one: 



Byx' U{t) = Ai) ^ \ X = MTV »=^h) (162) 
and hence: 

tu 2 + v! (1 + ut) = 0, (163) 

but as we supposed these attemps do not simplify our ode. 

Now we change the strategy and we are going to suppose that the ode has dimension of y. In this 
case we need to introduce the following dimensional constants in order to make eq. (158) verify the 
p.d.h. 

ay' = -bx 2 y 3 - cy 2 - a-, (164) 

x 

where [a] = x, [b] = y~ 2 x~ 2 and [c] = y , which bring us to the following c.v. 

= x, u(t) = ^\ =>• = t, y = -^—^J (165) 
therefore we obtain this new ode 

u ut - 2u 2 - t 5 - t 2 u = 0. (166) 

but as in the above tactic we have not advanced. This c.v. is precisely the suggested one by the theoretic 
methods. 

Following the Lie method we have to solve the pde: 

rjx + (Vv ~ £*) (~Ax 2 y 3 - By 2 - |) - £ y (-Ax 2 y z - By 2 - £)' - 

-i (J^y - 2Axy 3 ^j - V {-2By - ± - 3Ax 2 y 2 ^j = 0, (167) 

but this ode does not admit any symmetry. 

To end, we would like to show that the theoretical method does not work either. For this purpose 
we follow step by step the method beginning with a generic Abel ode of first order written as follows: 

y' = fsy 3 + hv 2 + hv + fo, (168) 

where in our case: 

h = -x 2 , / 2 = -l, /! = --, /o = 0. (169) 

x 

The c.v. suggested by the theoretical method is the following one 

x, u(t) = ^\^(x = t, y = -L^ , (170) 

(note that this c.v. is the same than the suggested one by D.A. (165)). In this way our ode is now 
rewritten as: 

u'u = h,2U 2 + h\u + ho, (171) 

where 

h 2 = -fi, h 1 = -f 2 , h = -f 3 , (172) 
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finding therefore that in our case we have: 

u 2 

u'u= — + u + t 2 , (173) 

which is an Abel of second order (this ode has no solution). To try to find a solution of this ode we 
make the following c.v. 

(r = t, s(r) = u(t)E) (t = r, u(t) = ^-J , (174) 

obtaining this new ode: 

ss' = F lS + Fo, (175) 

where 

E = exp(- J h 2 ), Fx = hiE, F = h E 2 , (176) 

hence 

ss' = - + 1, (177) 

r 

but unfortunately we do not know how to solve this apparently simple ode. 

As we have seen this ode seems very pathological since none of the followed tactics have helped us 
to obtain any solution. ■ 



5 Conclusions and discussion. 

We have seen how writing the odes in such a way that they verify the pdh i.e. introducing dimensional 
constants, we can obtain in a trivial way c.v. that bring us to obtain simpler ode than the original and 
therefore their integration is immediate. Furthermore, we have tried to show that these c.v. are not 
obtained as if by magic but that they correspond to invariant solutions or to particular solutions and 
therefore they are generated by the symmetries that admit the ode. 

Nevertheless, the D.A. has strong limitations. For example D.A. is unable (at least at this time 
we do not know how to do it) to solve the following simple linear ode 

y'=( x3 + l + 3 )y + ( 3x2 -^)> ( 178 ) 

for this reason one must not put all his confidence in this "tactic". This is one of the greater incon- 
venience that presents the proposed method. But as we have noticed in the introduction, we think 
that our pedestrian method continues having validity at least when one is studying ode derived from 
engineering problems or physical problems etc... where, as it is supposed, such odes must verify the pdh 
in such a way that for example the ode y' = x + 1, lacks of any sense (physical sense, since we cannot 
add a number to a physical quantity). 

Nevertheless, and in spite of the limitations that we have not avoided to show, we continued believing 
in the kindness (goodness) of the method and that it can be applied to obtain solutions (at least 
particular solutions and in concrete invariant solutions) to more complicated equations like the following 
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ones: 

y = ay, [a\=y x =^y = — , (179) 

y x z 

j/'" = 4' H = A~ 3 ^y = x 3/ \ (180) 

y"' = -ayy", [a] = y~ l x~ l =j- y = x~\ (181) 

y" = — +a^4, [a] = y~ x x =>• y = x, (182) 

= TTTZ , v [a]=y- 1 x^y = x, (183) 

y' (i + «y ) 

but these are questions that we will approach in a forthcoming paper. 
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into English. 

References 

[1] A.E. Ruak. Jour, of the Mitchell Socety. August (1935), 127-133. 

[2] G. Birhoff. Electrical Egineering. December (1948), 1185-1188. 
G. Birhoff. "Hydrodynamics" . Princeton (1950) 

[3] F.A. Morrison, Jr. Bull Mech. Eng. Edu. 8, (1969), 289-300. 

[4] M.J. Moran. Jour. Franklin Inst. 292, (1971), 423-432. 

[5] J. Hainzl. Jour. Franklin Inst. 292, (1971), 463-470. 

[6] T.Y. Na and A.G. Hansen. Jour. Franklin Inst. 292, (1971), 471-489. 

[7] G.I. Barenblatt and Y.B. Zel'dovich. (1972) 285-312 

[8] K. Kurt. "Dimensional Analysis and Group Theory in astrophysics" . Pergamon Press Oxford (1972). 

[9] R. Seshadri and T.Y. Na. "Group Invariance in Enginenering Boundary Value Problems". Springer- Verlang. NY. 
(1985) . 

[10] G.I. Bararenblatt. "Scaling, Self- Similarity and Intermediate Asymptotics" . CUP (1996) 

[11] M. Castans and C. Diaz. J. Non-Equilib. Thermodyn. 14, (1989) 315-329. 
M. Castans and C. Diaz. II Nuovo Cimento D 13, (1991), 1061-1069. 

[12] J. A. Belinchon and P. Davila. Class. Quantum Grav. 17,(2000), 3183-98. 
J. A. Belinchon . gr-qc/0404028. 

[13] M. Castans and P. Davila. GIAD. Lanzarote 1999. 

[14] E.J. Wilczynski. Trans. Am. Math Soc 23, (1900). 

[15] A. Cohen. "An Introduction to the Lie Theory of one-parameter groups" . Heath, NY. (1911). 

21 



[16] A.J.A. Morgan. Quart. Jour. Math. Oxford (2), 3 (1952) 250-59 

[17] L. P. Eisenhart. "Continuous Groups of Transformations". Dover 2003. 

[18] L.V. Ovsiannikov. "Group Analysis of differential Eqautions". Academic. Press, N.Y. (1982). 

[19] N.H. Ibra gimov. " Elementary Lie Group Analysis and Ordinary Differential Equations". Chicherter. John Wiley & 
Sons (1999). 

N.H. Ibragimov. "Introduction to Modern Group Analysis" . Ufa. (2000). 

[20] G.W.Bluman and J.D. Cole. "Similarity Methods for Differential Equations" . Springer- Verlang (1974). 

G.W.Bluman and S. Kumei. "Symmetries and Differential Eqautions". Springer- Verlang (1989). 

G.W.Bluman and S.C. Anco. "Symmetry and Integration Methods for Differential Equations". Springer- Verlang. 
(2002). 

[21] P.E. Hydon. Symmetry Methods for Differential Equations. CUP, (2000) 
[22] B.J. Cantwell. "Introduction to Symmetry Analysis". CUP. (2002). 

[23] G. Baumann. Symmetry Analysis of Differential Equations with Mathematica. Springer Telos. (2000). 
[24] H.Stephani. Differential Equations. Their Solution Using Symmetries. CUP (1989). 

[25] G. Emanuel. Solution of Ordinary Differential Equations by Continuos Groups. Chapma&Hall/CRC (2001). 

[26] J.M. Hill. Differential Equations and Group Methods forScintists and Engineers. CRC (1992). 

[27] P.J. Olver. Applications of Lie Groups to Differential Equations. 2nd ed. Springer- Verlang (1993). 

[28] E. Kamke. Differentialgleichungen: Losungsmethoden und Losungen, Akademische Verlagsgeselischaft Leipzig (1959). 

[29] G.M. Murphy. Ordinary Differential Equations and their Solutions. Princeton: Van Nostrand, (1960). 

[30] H.T. Davis Introduction to Nonlinear Differential and Integral Equation Dover N.Y. (1962). 

[31] D. Zwillinger. Handbook of Differential Equations. 2nd edition. AcademicPress, (1992). 

[32] A.D. Polyanin and V.F. Zaitsev. Handbook of Exact Solution for Ordinary Differential Equations. CRC Press Boca 
Raton FL (1995). 

[33] E.S. Cheb-Terrab and A.D. Roche. : European Journal of Applied Mathematics (2003), vol. 14, pp. 217-229. 
math.GM/0002059. 

E.S. Cheb-Terrab and A.D. Roche. Computer Physics Communications, math-ph/0001037 
E.S. Cheb-Terrab and T. Kolokolnikov. math-phys/0007023. 

[34] F. Schwarz. Studies in Applied Mathematics 100, (1998), 269-294. 

[35] V.M. Boyko. nlin.SI/0404020 

[36] M.K. Mak, H.W. Chan and T. Harko. Computers & Mathematics, 41 (2001) 1395-1401. 
M.K. Mak and T. Harko. Computers & Mathematics, 43 (2002) 91-94. 
T. Harko and M.K. Mak. Computers & Mathematics, 46 (2003) 849-853. 

[37] M.L. Green and P. A. Griffiths. Houston Journal of Mathematics 28, (2002) 329-351. 



22 



